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We investigate the properties of impenetrable bosons confined in a one-dimensional lattice at finite 
temperature in the presence of an additional incommensurate periodic potential. Relying on the 
exact Fermi-Bose mapping, we study the effects of temperature on the one-particle density matrix 
and related quantities such as the momentum distribution function and the natural orbitals. We 
found evidence of a finite temperature crossover related to the zero temperature superfluid-to-Bose- 
glass-transition that induces a derealization of the lowest natural orbitals. 



I. INTRODUCTION 

Recent developments in the field of ultracold atomic 
gases have opened up the possibility to investigating var- 
ious challenges from the physics of strongly correlated 
systems pQ. The remarkable experimental control over 
most of the system parameters, including dimensionality 
and interactions, makes it possible to use these setups as 
"quantum simulators" of fundamental Hamiltonian mod- 
els, otherwise inaccessible in more traditional condensed- 
matter realizations. For instance, the achievement of 
strong interactions and confinement to low-dimensional 
geometries led to the observation of the Tonks- Girardeau 
gas, with the addition of an optical lattice along the one- 
dimensional axis [2] and in the continuum [3]. 

A specially active line of research is the study of 
disorder-related phenomena [1]. One of the possibili- 
ties to create random potentials in a controlled way is to 
use laser beams generating speckle patterns [ST-flO]. An 
alternative technique is to superimpose two optical lat- 
tices with incommensurate frequencies, thus generating a 
quasiperiodic potential These realizations of ultra- 
cold atoms in disordered potentials enabled the observa- 
tion of Anderson localization in the regime of negligible 
interactions [THEI]. However, strictly speaking, Ander- 
son localization is a single-particle effect, and hence the 
quest turned to the understanding of the interplay be- 
tween disorder and interactions, a particularly intriguing 
aspect of the localization phenomenon on which recent 
experiments on bosonic systems [141 115] were performed. 
From the theoretical point of view, interacting bosons in 
a disordered potential are expected to exhibit a phase 
transition from a superfluid phase with extended single- 
particle states to a Bose-glass phase in which single- 
particle wavefunctions are Anderson localized as the in- 
teraction or the disorder strength is augmented [THl [17] • 
However, in spite of the theoretical efforts, the ground- 
state properties of disordered interacting bosons are still 
under debate [TM23] . and, for instance, the possibility of 
a reentrant superfluid induced by disorder [3] was exper- 
imentally excluded [15] . A finite temperature was pro- 
posed as a source of discrepancies [15]. The reason for 
this lack of agreement lies in the high complexity of the 
models analyzed such as the disordered Bose-Hubbard 



model, which can only be solved by means of sophisti- 
cated numerical techniques. In this paper we concentrate 
on the finite-temperature properties of interacting bosons 
confined on a one-dimensional quasiperiodic potential in 
the limit of infinite repulsion, for which an exact solu- 
tion exists |H] based on a mapping of the Hamiltonian 
to spinless fermions through the Jordan- Wigner trans- 
formation. We focus on the effect of the temperature 
on the behavior of off-diagonal one-particle correlations 
and related quantities such as the momentum distribu- 
tion and the natural orbitals and their occupations. We 
analyze the finite-temperature crossover in the vicinity of 
the transition from the superfluid to the Bose glass at a 
finite value of the strength of the secondary potential. We 
find that the lowest natural orbitals, which are localized 
on the Bose-glass phase at zero temperature, delocalize 
when the system enters into the quantum critical region. 



II. MODEL AND METHODS 

We consider hard-core bosons in a one-dimensional lat- 
tice subjected to a quasiperiodic potential, described by 
the following Hamiltonian: 
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where the boson destruction and creation operators bi 
and b\ satisfy the special commutation relations 



{b h b]} = l, 6f = (&f) 2 = 0, 



\h,b 



Mi 



[bib]} 



0. 



(2) 
(3) 



avoiding multiple occupation of a site. A determines the 
strength of the incommensurate potential, and a is an 
irrational number that characterizes the degree of incom- 
mensurability. Throughout this paper we treat the case 
of cr = (\/5 — l)/2, which is a common choice. 

The hard-core bosons Hamiltonian is exactly solvable 
by means of the Jordan- Wigner transformation, 
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that maps Eq. ([I]) onto an equivalent noninteracting 
fermionic Hamiltonian, 



A cos(27rjcr)n i , (5) 



where the c and satisfy the canonical anticommuta- 
tions relations. The boundary conditions must be han- 
dled with care since the last term of the sum over sites 
maps to b\b\ = — c' L c\e l7 ' N ^c* L ci, with N the opera- 
tor representing the total number of bosons. Thus, an 
even (odd) number of bosons induces antiperiodic (peri- 
odic) boundary conditions over the fermion model, which 
makes difficult its study in the grand-canonical ensemble. 
However, open boundary conditions are preserved by the 
mapping, and we consider this case in the present work. 

We are interested in the properties of the one-particle 
density matrix, defined as 
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At zero temperature (. . .) represents a iV-particle ground 
state expectation value. For this case, Rigol and Mu- 
ramatsu 25-27 found a very convenient way of repre- 
senting Eq. j6l, well suited for a numerical computa- 
tion. The one-particle density matrix can be written as 
Pij = Gij + 8ij(l — 2Ga), where is the one-particle 
Green's function and, according to Refs. [2"5Tf2"Tl can be 
represented in the form 



G ij = {b i b}) = det{Pl(i)P{j)} 



(7) 



The matrices P(i) are obtained by taking the first N rows 
from the matrix of the eigenvectors of Hp, then changing 
the sign of the first i — 1 rows and adding a column with 
vanishing elements, except for the i-th one, which is equal 
to 1. This provides an efficient way of computing pij 
since it involves the calculation of determinants of (N + 
1) x (N + 1) matrices instead of the more standard way 
that entails the evaluation of the Toplitz determinant of 
(L + 1) x (L + 1) matrices[2j [28]. One can take a step 
further and explicitly compute the product of matrices 
P^(i)P(j) for later reducing its determinant by using the 
property for block matrices 
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with Z the lxl matrix formed by the element N + 1,N + 
1 of the product P f (i)P(j). The one-particle Green's 
function can thus be expressed as 
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where ip m (i) are the wavefunctions of H. The N x N 
matrix A is 



A(i,j) = Q- 1 detQ, 



(10) 



where the entries of matrix Q are 
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(for j ^ i). The advantage of Eq. ([9| is that the nu- 
merical calculation of Gij for every pair of points can be 
optimized, and even parallelized, since Q(i,J + 1) can be 
easily updated from Q(i,j). Notice that Eqs. ([9|-([TT]) 
are the discrete version of the ones obtained in Ref. 1291 
for a one-dimensional system of impenetrable bosons in 
the continuum. 

At finite temperature T, Eq. ^ is obtained from 
a thermal average over a the grand-canonical ensemble 



described by a density operator p 



p -(H-»N)/k B T 



/z, 



where Z = tr e -( H -n N )/ k vT _ p^ s usu£ j ; th. e chemical po- 
tential p is determined implicitly by the condition that 
fixes the mean particle number N = pa. For the cal- 
culation of the nondiagonal elements we first employ the 
method proposed in Ref. [301 that makes use of the prop- 
erties of Slater determinants and leads to the following 
result for the one-particle density matrix: 

Pi j = |{dct[I + (I + A)0(j)URU t O^)] 

-detp + OyjUIttjtOCi)]}, (12) 

where U is the unitary matrix containing the eigenvector 
of H in its columns, 0(i) is diagonal with the first i — 1 
elements of the diagonal equal to — 1 and the others equal 
to 1, and the matrix A has only one non-zero element: 
Aij = 1. In addition, we defined R = e -(E-/ii)/fcBT 
with I the identity matrix and E = U^HU the diagonal 
matrix with the eigenvalues of H. 

However, the evaluation of these determinants can be 
unstable at low temperatures since the matrix R con- 
tains very small and very large elements in its diagonal. 
In order to obtain more stable numerical results, it is 
convenient to isolate R by rearranging the factors inside 
the first term of Eq. (fT2| as 



det [I + (I + A) 0(j)URU t O(i)] = det 0(j)0(i) 

x det [U t O(j)0(i)U + R + U t AU] , (13) 



and the same for the second term of Eq. (12 1. Moreover, 
we notice that the product U^AU can be written as an 
outer product, 



[TJ^AU], 



(14) 



Therefore, we can apply to the right-hand side of Eq. 
( 13 1 the matrix determinant lemma 



det (M + u ® v T ) = det M (l + v T M _1 u) 



(15) 



for any invertible square matrix M, and we finally obtain 
for the one-particle density matrix the result 

L 
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FIG. 1: Log-log plots of the decay of the one-particle density 
matrix with distance (in units of the lattice constant a) for 
ksT — 0.04t and different values of A in a half-filled systems 
with 1000 lattice sites. Thin solid lines are exponential fits. 
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where now the L x L matrix B is 



B(i,i) = (-1) 



; _.det(Q + R) 
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and (17) are the finite- 



and Q is the extension of the matrix defined in Eq 
to size L x L. In order to evaluate the one-particle den- 
sity matrix at finite temperature, one needs to calculate 
the L x L matrix Q + R, which can be updated from 
site to site, and its inverse and determinant for each pair 
of sites, for which there are fast and accurate numeri- 
cal procedures. Equations, 
temperature analogs of Eqs 



III. FINITE TEMPERATURE CROSSOVER 

At zero temperature, the phase diagram of the Hamil- 
tonian described by Eq. ([!]) has been extensively stud- 
ied [T8H2ol I3TI 132] , and it was found that three phases 
exist, namely, superfluid, Bose-glass, and incommensu- 
rate insulator, depending on the value of the chemical 
potential and A. One can relate this zero-temperature 
phase diagram to the single-particle band structure of the 
fermionic Hamiltonian of Eq. ^ . The energy spectrum 
of Hp is organized in a series of bands [33J EH] • Its wave 
functions are all localized for A > A c , with A c = 2t, 
and all extended for A < A c giving rise to the Bose-glass 
and the superfluid phases, respectively, when the chem- 
ical potential lies inside one of the bands. For values of 
the chemical potential inside a gap, the system is an in- 
commensurate insulator for any value of A. By fixing 
the chemical potential in the middle of a band, one can 
study the transition from the Bose glass to the superfluid 
by varying A. We explicitly analyze this case by setting 
fx = in all the calculations which fixes the mean number 
of particles to half filling. 



FIG. 2: Correlation length as a function of temperature and 
A. (a) f(T) for A = t (in logarithmic scale). Solid thin lines 
show power law fits, (b) for A = 3.5t in log-log scale. A 
power-law behavior is observed for ksT > 0.25i. (c) £(A) for 
ksT — 0.05t. (d) Same as (b) in log-linear scale. 



At zero temperature the one-particle density matrix 
exhibits a power law decay pij cx |xj — Xj^ 11 for A < A c , 
characteristic of a quasi-long-range ordered superfluid 
phase. In the absence of the quasiperiodic potential, 
the exponent is 77 = 1/2 35J. As the incommensu- 
rate secondary potential is turned on r\ becomes smaller, 
wwhich is accompanied by a quick decrease of the su- 
perfluid fraction [32]. When A exceeds the critical point 
A c , the one-particle density matrix decay turns into an 
exponential law ~ e -\ x i- x i\/£ j where £ is the corre- 
lation length, signaling the transition to the disordered 
Bose-glass phase [T7]. In one-dimensional (ID) systems 
any finite temperature destroys quasi-long-range order, 
inducing an exponential decay of the correlations char- 
acterized by a temperature- (and A in the problem un- 
der study) dependent correlation length. For hard-core 
bosons in the absence of the quasiperiodic potential this 
has been verified earlier [30] , where the power-law behav- 
ior of the correlation length with increasing temperature 
£ ~ 1/T was found for low temperatures k^T < t. Here 
we extend these calculations to finite A, investigate the 
finite-temperature crossovers associated with the quan- 
tum critical point at A = A c [36] . and characterize the 
properties of the different regions in the T — A phase 
diagram. 

Figure [l] shows the decay of the one-particle density 
matrix for a half- filled system of 1000 lattice sites at a 
finite temperature k^T — O.OAt and different values of 
A. The decay is exponential, characterized by a corre- 
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FIG. 3: Finite T phase diagram of ID hard core bosons in a 
quasiperiodic optical lattice at half-filling in the A — T space 
in the vicinity of the quantum critical point at A = A c = 2t. 
Quasi-long-range order is present at T — for A < A c as 
a superfluid (SF). The ground state for A > A c is a Bose 
glass (BG). Regions I, II, and III are characterized by the 
temperature dependence of the correlation length. We sketch 
the phase boundary between regions II and III. 



FIG. 4: Estimated exponents of the power-law behavior of the 
correlation length, £ ~ T _M , as a function of A. For A < A c 
we show the low-temperature exponent /i_ in brown (solid) 
and the higher-temperature exponent [i+ in blue (dashed). In 
orange (dotted) we represent the high-temperature exponent 
for A > A c . Notice that at A = A c both low- and high- 
temperature exponents tend to converge to a common value. 



lation length £ whose temperature and A dependence is 
shown in Fig. [2] The study of the behavior of £ for dif- 
ferent regimes allows us to establish the crossover phase 
diagram shown in Fig. [3] For values of A ^ A c , we find 
that the correlation length diverges at low temperatures 
as a power law £ ~ [Fig.[2^a)], with an exponent /x_ 

that depends on A. This is marked as region II in Fig. [3] 
and can be identified with a thermally disordered phase 
in which quasi-long-range superfluid order is destroyed 
by thermal fluctuations. We also observe a crossover to 
a power law behavior with a different exponent /i+ for 
temperatures above a crossover temperature XI (A), (re- 
gion III in Fig. [3]). The obtained data allows us to make a 
sketch of the phase border between regions II and III. On 
the localized side of the spectrum the correlation length 
behaves as a power law with increasing temperature for T 
larger than the right branch of the crossover temperature 
T+(A) [Fig[2^b)], but when lowering the temperature, 
it saturates to the zero-temperature finite value follow- 
ing the exponential behavior £ ~ e~ aT for T < T + (A) 
[Fig|2td)]. This is represented in Fig. [3] as region I, and 
can be identified as a quantum disordered phase, since 
its finite-temperature properties are dictated by the lo- 
calized Bose-glass ground state. In particular, as we see 
in the next section, natural orbitals are still localized in 
this regime. Region III, in between I and II, is essentially 
the quantum critical region in which both types of fluc- 
tuations, quantum and thermal are important. A more 
precise determination of the phase borders, in particular 
between phase II and III, will be done elsewhere 37J. 

Figure [4] shows the nonuniversal behavior of the ex- 
ponents n± as a function of A for the whole range 



of values of A and the two temperature regimes when 
A < A c . For A < 1.5t the exponent /z_ is almost con- 
stant retaining the zero A value /i_ = 1 whereas for 
A > 1.5f it decays abruptly to its limit value at A = A c 
(for A > A c the low-temperature behavior changes to 
an exponential form). On the other hand, the high- 
temperature exponent /i + shows a smoother decreasing 
behavior with A, that continues to the A > A c region. 
At A = A c , we still have two regimes with the same 
exponent /i + ~ /i_ ~ 1/2 that differ in the prefactor. 
Even though the low-temperature region is very small, 
we should expect a single high-temperature regime at 
the critical point. We attribute the survival of the low- 
temperature one to a finite size effect. 

Figure [2^c) shows the dependence of £ on A for k^,T = 
0.05i. For low A the correlation length converges to the 
zero disorder value. As A is raised the correlation length 
smoothly switches to a power law at a certain value of the 
secondary lattice amplitude that moves to smaller values 
for increasing temperature. For even higher amplitudes 
of the incommensurate potential we find a new crossover, 
and the behavior departs from a power law. 

The behavior of the one-particle density matrix for 
A < 1 and fillings such that we are far from the band 
edges can be understood in terms of the Luttinger liquid 
description of the ID interacting bosonic fluid, using a 
low-energy hydrodynamic approach [THl HOI I58H4T3] . In 
particular, the order parameter correlation function in 
the thermodynamic limit (and in the continuum) is given 
by gD] 

\smh(7r|a; — x'\/ Lt) J 
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FIG. 5: Momentum distribution for A = 1.9t (main plot) 
and A = 2.5t (inset) at different temperatures for a half- filled 
system of 800 lattice sites. At A = 2.5t curves for k^T — 
and ksT = O.Olt are identical. 



FIG. 6: Condensate fraction as a function of A for different 
temperatures. At A = A c it reaches a common residual value 
signaling the transition to the Bose glass. 



where Lt ~ 1 /T is the thermal correlation length and K 
is the Luttinger liquid parameter, which takes the value 
K = 1 in the absence of the secondary lattice, but is 
renormalized at finite A. From this result one observes 
that correlations cross over from algebraic decay for |a; — 
x'\ << Lt with exponent 1/2K to an exponential form 
for |ar — a;'] >> Lt- The latter case results in a behavior 
£ r~> 1/T in agreement with Ref. [30] and the results of Fig. 
[I] The fact that the exponent /Lt departs from the value 
\i = 1 for A > t signals the breakdown of the Luttinger 
liquid description as the system approaches the quantum 
critical region. 

Let us finally analyze the momentum distribution func- 
tion, 

ij 

where k denotes the momentum. In Fig. [5] we show the 
behavior of the momentum distribution for two represen- 
tative values of A in a system of 800 lattice sites. A finite 
temperature reinforces the effect of the pseudodisorder on 
the correlations. At A = t, temperature washes out the 
two fundamental features of the momentum distribution: 
the central peak (associated with the quasi-long-range 
correlations at T = 0) and the secondary peaks (related 
to the secondary incommensurate lattice). In the ab- 
sence of the secondary lattice, at k-gT = O.Oli the height 
of the central momentum peak was found to decrease to 
2/3 of the one at zero temperature [30]. A finite value of 
A = 1.9t at the same temperature reduces even further 
its height to around 1/3 of the central one, though the 
secondary peaks are still visible. At k^T = 0.5t both 
features essentially disappeared. 



IV. DELOCALIZATION OF THE NATURAL 
ORBITALS 

Since the single-particle wavefunctions of the Hamilto- 
nian ^ are all either localized for A > A c or delocalized 
for A < A c , that is, there are no mobility edges, one ex- 
pects activation effects to be suppressed and temperature 
to play no role in the physics of derealization. However, 
the Jordan- Wigner transformation leads to a highly com- 
plicated nonlocal form of the bosonic operators, which 
keeps track of the strong interactions, and might account 
for nontrivial effects in the vicinity of a quantum critical 
point, hidden in the apparent simplicity of the fermionic 
Hamiltonian. We analyze here these effects by studying 
the localization properties of the natural orbitals defined 
as the eigenvectors of the one-particle density matrix: 

JV 

5>^ = A„C- (20) 

3=1 

Here the corresponding eigenvalues {A"} represent the 
occupations. The natural orbitals can be considered as 
effective one-particle states [41] in systems with interac- 
tions. For noninteracting systems the natural orbitals 
reduce to the exact independent-particle cigenfunctions 
of the Hamiltonian, and at finite temperature their oc- 
cupancies are either the Bose or Fermi factor depend- 
ing on the character of the statistics. However, when 
interactions are present as in Eq. (JlJ, and unlike the 
eigenstates of the Hamiltonian the natural orbitals 
depend on temperature |30j . In dilute high-dimensional 
gases exhibiting Bose-Einstein condensation (BEC), the 
lowest natural orbital is usually regarded as the conden- 
sate wave function, and the BEC order parameter [42] . 
In systems with disorder, the phenomenon of Anderson 
localization in dilute gases is usually studied by consider- 
ing the localization of the condensate wave function [33] . 
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Here we extend this analysis to higher natural orbitals 
since they turn out to be occupied when the system ap- 
proaches the Bose-glass phase. 

Let us investigate first the occupations of the natu- 
ral orbitals. The superfluid phase at zero temperature 
is associated with the formation of a quasi-condensate, 
that is, a large occupation of the lowest natural orbitals 
can occur. For instance, for hard-core bosons in the ab- 
sence of the secondary lattice, Ai scales as ~ VN [2"6] . 
In Fig. [6] we observe the occupation of the first orbital 
for different temperatures along with the T — results 
as a function of the incommensurate potential amplitude. 
For all temperatures the condensate occupation decreases 
with increasing A, and reaches a common residual value 
at A = A c the value at which the single-particle eigen- 
functions localize and the system enters the Bose-glass. 
Another signature of the Bose-glass phase at zero temper- 
ature is the opening of a discontinuity in the occupation 
of the orbitals for v = N, resembling the Fermi-Dirac 
distribution, as can be observed in Fig. [7J in which we 
show the distribution of occupations of the orbitals versus 
the orbital number for fixed A at different temperatures. 
However, the analogy with the Fermi-Dirac distribution 
function is not complete, since the occupations of natural 
orbitals below a certain threshold (of around v = N/2) 
become larger than 1 and slowly increase when descend- 
ing in the orbital number. The discontinuity at v = N 
is erased by thermal effects even at very low temper- 
atures. However, at higher temperatures (k^T ~ 0.5i 
for A = 2.5t) two discontinuities develop symmetrically 
placed around v = N. At very high temperatures we 
approach the infinite-temperature regime, in which all 
occupations are equal to 1/2. 

At zero temperature, the superfluid-to-Bose-glass tran- 
sition is marked by the localization of the lowest orbitals 
at A = A c . The insets of Figs.^a) and^b) show the ef- 
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FIG. 7: Occupations of the natural orbitals at A = 2.5t and 
different temperatures for a half-filled system of 800 lattice 
sites. The abrupt jump observed at v — 400 at zero temper- 
ature is rounded as temperature rises. For very high temper- 
atures all the orbitals are equally occupied. 



feet of the incommensurate potential on the lowest orbital 
below and above the transition. For A < A c the states 
are extended, whereas for A > A c they exponentially lo- 
calize, which mimics the behavior of the single-particle 
states of the fermionic Hamiltonian of Eq. Never- 
theless, as one approaches the point A = A c the range of 
significantly occupied high-order orbitals increases, im- 
plying that the dominance of the lowest ones decreases. 
Therefore, the characterization of the system by a single 
macroscopic wave function <f>x can be put into question 
with increasing A, and one is led to analyze the whole set 
of significantly occupied orbitals. Given a single-particle 
state \4>) we can quantify its degree of localization by the 
inverse participation ratio (IPR) in real space, defined as 

L 

ipr(<^]T|(^)i 4 (21) 

where \ j) is the state at site j in the real-space basis. The 
IPR is a measure of the inverse number of sites occupied 
by the state as can be understood by regarding two lim- 
iting cases. If the state is localized on a single site i, we 
immediately find from (<f>\j) = #y that IPR = 1. If, on 
the other hand, the state is distributed equally over L 
sites, we have (<f>\j) = L^ 1 / 2 and thus IPR = 1/L. 

Figures |8ja) and |8jb) show the IPR for the whole 
set of natural orbitals together with their occupation for 
two values of A, below and above the transition at zero 
temperature. We observe that for A < A c , where all 
the single-particle fermionic states are extended, there is 
a group of localized natural orbitals between v = 300 
and v = 600 for a half filled system of 800 lattice sites. 
For A > A c , in which case all single-particle fermionic 
wave functions are localized, some natural orbitals are 
extended between v = 50 and v = 150. 

The effect of a finite temperature on the natural or- 
bitals is shown in Figs. |8^c) and J8jd) where the profile 
of the first natural orbitals for A = 3t at two different 
temperatures is plotted in the insets. As the tempera- 
ture is increased from zero to a small but finite value, 
the lowest natural orbitals are still exponentially local- 
ized. This is in agreement with the fact that in this 
finite-temperature region, the physics is dominated by 
the quantum fluctuations of the Bose-glass ground state. 
When the temperature exceeds some value T4 = T^(A), 
of about Td ~ 0.05t for A = 3t the lowest natural orbitals 
are not localized anymore; instead they spread over all 
the lattice sites, even for A > A c . This suggests a sort 
of competition between the effects of the incommensu- 
rate potential that acts like a random potential on sites 
because of the irrational value of a and tends to localize 
the natural orbitals and thermal effects traceable solely 
trough correlations of nonlocal operators for which the 
system is fully interacting, that aims to delocalize them. 
Nevertheless, as in the zero-temperature case, not all the 
orbitals exhibit the same behavior, not even all the occu- 
pied ones. Figure [8] also shows the occupations and the 
IPR of all the natural orbitals, and we observe the same 
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FIG. 8: Occupation and IPR of the orbitals at (a) k B T = and A = t, (b) k B T = and A = 3t, (c) k B T = O.Oli and A = 3t, 
(d) k B T = O.lt and A = 3t for a half-filled system of 800 lattice sites. In the insets we show the profile of the first orbital. 
Note the log scale in (b) and (c). 



situation as at zero temperature: The derealization oc- 
curs for orbitals such that v < 300, whereas the group of 
orbitals around and above the step in the occupations at 
v = N stay localized. Finally, we show in Fig. [9] the IPR 
of the first orbital as a function of temperature for dif- 
ferent values of A. For the two values of A > A c shown, 
the IPR reduces to a residual value when the tempera- 
ture exceeds T^(A). Remarkably, the temperature of 
derealization of the lowest orbitals is very close to the 
crossover temperature T + at which the correlation length 
changes its behavior as a function of temperature. 



V. SUMMARY 

In conclusion, we have studied finite temperature prop- 
erties of hard-core bosons in a ID incommensurate lattice 
potential. Based on the Jordan- Wigner transformation 
and by using properties of Slater determinants, we com- 
puted the one-particle density matrix and related quan- 
tities such as the momentum distribution and natural 



0.5 
0.4 
0.3 
0.2 
0.1 






0.1 



0.2 

k B Tlt 



0.3 



0.4 



FIG. 9: Inverse participation ratio of the first orbital as a 
function of temperature for different values of the secondary 
lattice amplitude. (Inset) A dependence of the IPR for the 
first orbital at T = 0. 
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orbitals for a half-filled system. We obtained a finite- 
temperature crossover phase diagram in which different 
regions are classified according to the behavior of the 
correlation length with temperature. For values of the 
incommensurate potential amplitude below the critical 
value A c we find that the correlation length diverges at 
low temperatures as a power law with a A-dependent ex- 
ponent. As the system enters into the quantum critical 
region, the correlation-length crosses over to a power-law 
form with a different exponent. In contrast, for A > A c 
there is a crossover from a low-temperature exponential 
behavior of the correlation length, characteristic of the 
quantum fluctuations in the Bose-glass ground state, to 
the power-law form when the systems enters into the crit- 



ical region on the localized side of the spectrum. We 
also found evidence that this crossover is characterized 
by an important change in the localization properties of 
the natural orbitals, in particular a significative fraction 
of the natural orbitals delocalize when the temperature 
exceeds the crossover temperature scale. 
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